Dynamics of continued fractions and kneading sequences of 

unimodal maps 

C. Bonanno* C. Carminati^ S. Isola * G. Tiozzo§ 

June 1, 2012 



Abstract 

In this paper we construct a correspondence between the parameter spaces of two families of 
one-dimensional dynamical systems, the a-continued fraction transformations T a and unimodal 
maps. This correspondence identifies bifurcation parameters in the two families, and allows one 
to transfer topological and metric properties from one setting to the other. As an application, 
we recover results about the real slice of the Mandelbrot set, and the set of univoque numbers. 



1 Introduction 

The goal of this paper is to discuss an unexpected connection between the parameter spaces of two 
families of one-dimensional dynamical systems, and establish an explicit correspondence between 
the bifurcation parameters for these families. 

The family (T a ) a e(o,i] °f a-continued fraction transformations, defined in |Na] . is a family of 
discontinuous interval maps, which generalize the well-known Gauss map. For each a G (0, 1], the 
map T a from the interval [a — 1, a] to itself is defined as T a (0) = and, for i/O, 



where c a . 



is a positive integer. For every a £ [0, 1] at most two of the infinitely 
many branches of T a fail to be surjective, hence T a has infinite topological entropy. On the other 
hand, for all a > the map T a admits a unique invariant probability measure \x a which is ab- 
solutely continuous with respect to Lebesgue measure, and one can consider the metric entropy 
h = h(T a , Ha) with respect to this measure; h can be computed by Rohlin's formula and is finite 
(see [EMJ). 
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Several authors have studied the variation of h as a function of the parameter a. Numerical 
evidence that the entropy h is continuous but non-monotone was first produced in [LM]. Subse- 
quently, |NN] found out that the entropy is monotone over intervals in parameter space for which 
the orbits of the two endpoints collide after a finite number of steps (see Equation ^ on page [6]) . 
This analysis is completed in |CT| , where intervals of parameters for which such relation holds are 
classified, and it is proven that the union of all such intervals has full measure. The complementary 
set, denoted by £ , is the set of parameters across which the combinatorics of T a changes, hence it 
will be called the bifurcation set. 

The second object we consider is the family of unimodal maps, i.e. smooth maps of the interval 
with only one critical point, the most famous example being the logistic family. To any such map 
one can associate a kneading invariant |MT| which encodes the dynamics of the critical point and 
determines the combinatorial type of the map. Using an appropriate coding (see jlPJ, [Is] ) the set 
of all kneading invariants which arise from unimodal maps can be represented as the points of a 
set A defined in terms of the tent map T(x) := min{2x, 2(1 — x)} 

A := { x £ [0, 1] : T k (x) < x Vk £ N} (1) 

The set A is uncountable, totally disconnected and closed; the topology of A reflects the variation 
in the dynamics of the corresponding maps as the parameter varies: for instance families of isolated 



points in A correspond to period doubling cascades (see Section 4.2). This set also parametrizes the 
real slice of the boundary of the Mandelbrot set Ai (see Section 5.1), because for each admissible 
kneading sequence there exists exactly one real parameter on the boundary with that kneading 
sequence. 

The key result of this paper (Section [4]) is the following correspondence between the bifurcation 
parameters of the two families: 

Theorem 1.1. The sets A\{0} and£ are homeomorphic. More precisely, the map ip : [0, 1] — > [J, 1] 

given by 

' — *+ip(x) = 0.11... 100... 011... 1... 



1 

a 2 + 



ai H j ai a2 <23 



1 

03 H 



is an orientation-reversing homeomorphism which takes £ onto A\ {0}. 

As a consequence, intervals in the parameter space of a-continued fractions where a matching 
between the orbits of the endpoints occurs are in one-to-one correspondence with real hyperbolic 
components of the Mandelbrot set, i.e. intervals in the parameter space of real quadratic polynomi- 
als where the orbit of the critical point is attracted to a periodic cycle. In terms of entropy, intervals 
over which the metric entropy of a-continued fractions is monotone (and conjecturally smooth) are 
mapped to parameter intervals in the space of quadratic polynomials where the topological entropy 
is constant (see figure [l]). For instance, the matching interval ([0; 3], [0; 2, 1]), identified in [LMJ and 
[NNJ, corresponds to the "airplane component" of period 3 in the Mandelbrot set. 

This dictionary illuminates many connections between seemingly unrelated objects in real dy- 
namics, complex dynamics, and arithmetic, which we will explore in the rest of the paper. For 
example, the following properties of A follow immediately by using the combinatorial tools devel- 
oped in [CT] : 
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Figure 1: Correspondence be- 
tween the parameter space of 
a-continued fraction transforma- 
tions and the Mandelbrot set. 
On the top: the entropy of a- 
c.f. as a function of a, from 
[LMJ; colored strips correspond 
to matching intervals. At the 
bottom: a section of the Mandel- 
brot set along the real line, with 
external rays landing on the real 
axis. Matching intervals on the 
top figure correspond to hyper- 
bolic components on the bottom. 



(i) the set A can be constructed via a bisection algorithm (Section 4.1) 

(ii) the sequences of isolated points in £ observed experimentally in ( |CMPT] . Section 4.2) are 



the images of the well-known period doubling cascades for unimodal maps (Section 4.2) 
(iii) the derived set A' is a Cantor set (Section |4.3[) 



(iv) the Hausdorff dimension of A is 1 (Section 4.4) 



From the above properties of A we derive some consequences on corresponding sets. For exam- 
ple, (iv) implies that the set of external rays of the Mandelbrot set which land on the real axis has 
full Hausdorff dimension (Section |5.1[ ), a result first obtained in [Za]. 

Moreover, in Section |5.2| we give an arithmetic interpretation of A. Namely, aperiodic elements 
of A correspond to binary expansions of univoque numbers, i.e. the numbers q G (1,2) such that 1 
admits a unique representation in base q. Univoque numbers have been studied by Erdos, Horvath 



and Joo |EHJj and many others (see Section 5.2 for references), and once again the translation of 
statements (i)-(iv) immediately yields several results previously obtained by different authors. 



Moreover, our new characterization of the set £ (see Lemma 3.3) shows that it is essentially the 



same object appearing in a paper of Cassaigne as the spectrum of recurrence quotients for cutting 
sequences of geodesies on the torus ( |Caj . Theorem 1.1). Indeed, our results on £ answer some 
questions raised in |Caj . such as the computation of Hausdorff dimension. 

Finally, we would like to emphasize that the correspondence of Theorem |1.1| opens up many 
questions, and it is especially natural to ask to what extent results in the well-developed theory 
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of unimodal maps can be translated into the continued fraction setting. For instance, it would 
be interesting to identify an analogue of renormalization for the a-continued fractions, and to 
further explore the relation between the entropy of that family ([LM|, [CM PTj . |Tij . |KSS] ) and 
the topological entropy of unimodal maps ( |MT| and |Doj among others). 
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2 Preliminaries 

Let T, F, G denote the tent map, the Farey map and the Gauss map of [0, 1], given b}Q 

2x 



T(x) 



if < x < \ 



2(1 - x) if \ < x < 1 



F(x) := < 



if < x < \ 

1 — x z 



1 — X 



if | < x < 1 



and G(0) := 0, G{x) := {^} , x ^ 0. The action of F and T can be nicely illustrated with 
different symbolic codings of numbers. Given x £ [0, 1] we can expand it in (at least) two ways: 
using a continued fraction expansion, i.e. 



x 



= [0;ai,a 2 ,a 3 , . . .] , ajGN 



ai + 



Cl2 + 
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and a binary expansion, i.e. 



Y,h2- l = o.6ib 2 ... , he {0,1} 



i>l 



The action of T on binary expansions is as follows: for uj S {0, 1} N , 

T(0.0w) = 0.w , T(0.1a;) = 0.w 



(2) 



where & = 0)10)2 . . . and 6 = 1, 1 = 0. The actions of F and G are given by F([0; 01, 02, 03, • • • ]) = 
[0;ai - l,a 2 ,a 3 ,...] if a\ > 1, while F([0; 1, a 2 , a 3 , ...]) = [0; a 2 , a 3 , ...], and G([0; 01, a 2 , a 3 , . . . ]) = 
[0; a 2 , a 3 , . . . ]. As a matter of fact G can be obtained by F by inducing on the interval I\ = [1/2, 1], 
i.e. 

G(x) = F^(x) , x^0 (3) 



1 Here [a;J and {x} denote the integer and the fractional part of x, respectively, so that x — [x\ + {x}. 
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Now, given x = [0; a±, G62, a>3, • • • ]> one may ask what is the number obtained by interpreting the 
partial quotients a>i as the lengths of successive blocks in the dyadic expansion of a real number in 
[0, 1]; this defines Minkowski's question mark function ? : [0, 1] — >• [0, 1] 



= 2 -(«i+-+«fc- 1 ) = 0. 00 ... 11 ... 1 00 ... 



k>l 



ai— 1 



a 2 



«3 



The questionmark function ?(x) has the following properties (see [SaJ): 



(4) 




(a) Tent map 



(b) Minkowski map 



(c) Farey map 



• it is strictly increasing from to 1 and Holder continuous of exponent f3 = — lo ^ +1 ; 

• x is rational iff ?(x) is of the form k/2 s , with k and s integers; 

• x is a quadratic irrational iff ?(x) is a (non-dyadic) rational; 

• ?(x) is a singular function: its derivative vanishes Lebesgue-almost everywhere; 

• it satisfies the functional equation ?(x)+?(l — x) = 1. 

We shall see later that the question mark function ? plays a key role in the main result of this 



paper (Theorem 1.1) because it conjugates the actions of F and T. For a more general analysis of 
the role played by ? in a dynamical setting see |BI] . 



3 The bifurcation sets 

3.1 The bifurcation set for continued fraction transformations 

Let r G (0, 1)DQ be a rational number, and let r = [0; a±, . . . , a n ] be its continued fraction expansion, 
with a n > 2. The quadratic interval associated to r will be the open interval I r with endpoints 



[0;ai,. . . ,a n _i,a n ] and [0; at, . . . , a„_i, a n - 1, 1] 

where [0; a\, . . . , a n ] denotes the real number with periodic continued fraction expansion of period 
(ai, . . . , a n ) . The number r is called the pseudocenter oiI r . We also define the degenerate quadratic 
interval I\ := (g, 1], where 5 := [0; 1] is the golden number. 
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Definition 3.1. The bifurcation set for continued fractions is 

£:=[0,1}\ |J I r (5) 

re(o,i]nQ 

Quadratic intervals which are not properly contained in any larger quadratic interval are called 
maximal quadratic intervals. By ( [CTj . Section 2.2), if two quadratic intervals are not disjoint, 
they are both contained in a maximal quadratic interval; thus, the connected components of the 
complement of £ are precisely the maximal intervals. For all parameters belonging to a maxi- 
mal quadratic interval, the ct-continued fraction transformation T a satisfies a matching condition 
between the orbits of the endpoints, with fixed combinatorics: 

Theorem 3.2 ([CTJ, Theorem 3.1). Let I r be a maximal quadratic interval, and r = [0; a±, . . . , a n ] 
with n even. Let 

N = ^2 ai M = ^2 ai 

i even i odd 

Then for all a G I r , 

T^\a)=T^\a-l) Va € I r (6) 

As a consequence, by |NN| . the metric entropy h(T a ) is locally monotone on the complement 
of £. The set £ can be given the following new characterization: 

Lemma 3.3. Let x G [0, 1]. The following are equivalent: 

(a) x £ £ 

(b) G k (x) > x, V/c G N 

(c) F k (x) > x, V/c G N 

(d) F k o ^i(x) < ^i(x) V/c G N, where ^\{x) := 1/(1 + x) is an inverse branch of F. 

Proof. Let x = [0; ai, a%, . . . ] be the continued fraction expansion of x. 

(b) => (c). Fix k > 1, so F k (x) = [0; /i, a n +i, a n+ 2 5 • • ■ ] for some re and 1 < h < a n . Hence 

F k (x) > [0;a n ,a n+1 ,...] =G n ~ 1 (x) >x 

(c) (d). Let J C N be the subsequence given by the partial quotients sums, i.e. J = 
{Yli=i a i}n>i, and note that y := ^i(x) = [0; 1, ai, a 2 , . . . ]. If fc ^ J then F fc (y) < i < y. 
Otherwise, F k (y) = [0; 1, a n+ i, a n+ 2, • • • ] = x I / i(i ?fc (x)) for some n > 1 so that 

x = [0; 01,02, a 3 , ...} < [0; a n+ i,o n+2 , . . .] = F k (x) F k (y) < y 

(d) (6). Fix k > 1, and let s := Oj, so that F s (x) = G fc (x). Now, 

F s (^i(x)) = [0;l,a fc+ i,a fc+2 ,...] < [0; 1, a u a 2 , . . . ] = ^i(x) 
therefore [0; a^+i, . . . ] = G k (x) > x = [0; ai, . . . ]. 

(a) =^ (c). Let x G £ be the endpoint of a maximal interval, hence x = [0; ai, . . . , a n ]. Once you fix 
k > 0, then G fc (x) = [0; aj+i, • ■ • , a n , ai, . . . , a/] with i = k mod re. Now, by [CT], Proposit ion 4.o 

G fc (x) = [0;aj + i, . . . ,a n ,ai, . . . , ai] > [0;ai, . . . ,a z ,a;+i, ... ,a n ] = x 
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i.e. (6), and hence (c). Since £ has empty interior ( |CT| . Theorem 1.2), every point in £ is a 
limit point of a sequence of endpoints of maximal intervals, so inequality (c) extends to all of £ by 
continuity of F. 

(b) => (a). If x ^ £, then x belongs to a maximal interval I a = (a - , a + ). Let a = [0; A~] = [0; j4 + ] 
be the two continued fraction expansions of the rational number a, in such a way that a~ = [0; A~] 
and a + = [0; ^4 + ]- Here A~ is a finite string of integers of odd length i, and A + is a finite string of 
integers of even length m. The sequence 

a = [0; A+] < [0; A+A+] < • • • < [0; (A+) n ] < [0; (^+) n+1 ] < . . . 

tends to [0; ^4+] = a + , hence for any x G [a, a + ) there exists n > such that 

[0;(^ + ) n+1 ] < x < [0;(^+) n+2 ] 

therefore y := G m (x) < [0; < x. If instead x £ (a - , a], then using the fact that / is odd 

= G l (a) < G l (x) < G\oT) = oT 

so G l (x) < a~ < x. □ 



3.2 Universal encoding for unimodal maps 

We now introduce a set A which encodes the topological dynamics of unimodal maps in a universal 
way. Several similar approaches are possible ( jMTj , |dMvS] among others) ; we follow jlP] , [ii] . We 
recall that a smooth map / : [0, 1] — > [0, 1] is called unimodal if it has exactly one critical point 
< co < 1, and we are going to assume /(0) = /(l) = 0. If x G [0, 1] never maps to cq under /, let 
us define the itinerary of x G [0, 1] to be the sequence 



■/ ^ f iff J W < co 

i{x)=sis 2 ... with Si = < - f H-U w 

[ 1 if p 1 (x) > c 



If x eventually maps to co, let us define i(x^) := lim^^i i(y) where the limit is taken over non- 
precritical points. The kneading sequence of a unimodal map is the binary sequence K := i(c[), 
where c\ = /(co) is the critical value. 

Kneading sequences encode the combinatorics of the orbits and determine the topological en- 
tropy. A nice way to decide whether a given binary sequence s G {0, 1} N is the kneading sequence 
of some unimodal map is the following: we first associate to s = s\S2 ■ ■ ■ G {0, 1} N the number 
t(s) G [0, 1] defined as 

k 

T = 0.t 1 t 2 --- = J2 t k2~ k , tfc = Si (mod 2) (7) 

k>l i=\ 

One readily checks that the following diagram is commutative: 

[0,1] — U {o,i} N [0,1] 



T (8) 



[0,1] — {0,1} N — !-)■ [0,1] 
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where T is the tent map (J2J) and a : {0, 1} N — > {0, 1} N is the shift map. Since the critical value 
is the highest value in the image and the map r o i is increasing ([Is], Lemma 1.1), one has the 
characterization 

Proposition 3.4 ( |IPj . [Is] ) . vl binary sequence s £ {0, 1} N is the kneading sequence of a continuous 
unimodal interval map if and only if t := t(s) satisfies 

T k (r) < t for all k > 

One is thus led to study the set 

A : = { x G [0, 1] : T k (x) < x, V7c G N} (9) 

By considering T{x) < x one immediately realizes that A \ {0} C [2/3, 1]. Note moreover that, if 
the map is C l , the extremal situations in which T fc (r) = r for some k > correspond to periodic 
attractors (and thus periodic kneading sequence) of the corresponding smooth unimodal maps (see 
[dMvS], Chapter IV). 

Note that all kneading sequences which arise from unimodal maps can be actually realized by 
real quadratic polynomials f c (z) = z 2 + c, c S [—2, j]. By density of hyperbolicity ( |GSj . |Ly| ), 
aperiodic kneading sequences are realized by exactly one map f c . On the other hand, the set of 
parameters which realize a given admissible periodic kneading sequence is an interval in parameter 
space, and among those maps there is exactly one f c which has a parabolic orbit. In conclusion, 
A parametrizes the set of topologically unstable parameters, hence it will be called the binary 
bifurcation set. 



4 From continued fractions to kneading sequences 



We are now ready to prove Theorem 1.1 namely that the map (p : [0, 1] — > 1] given by x = 
[0; ai,Q2, as, ...} h-> ip{x) = 0. 11 ... 1 00 ... 11 ... 1 . . . is an orientation-reversing homeomorphism 



«i 



a 2 



II:', 



which takes £ onto A \ {0}. 

Proof of theorem ] 1.1\ The key step is that Minkowski's question mark function conjugates the 
Farey and tent maps, i.e. 

l(F{x)) = T(?{x)) Vx 6 [0, 1] (10) 

Note that ip =? o hence (p is a homeomorphism of [0,1] onto the image of i.e. [^,1]- 
Moreover, by Lemma 3.3 -(d), 



#i(£) = {x e (0,1] : F k (x) < x} 



and by ( 10 ) 



?(*i(£)) = {x£ (0, 1] : T k (x) < x} = A \ {0}. 



□ 



In the following subsections we will investigate a few consequences of such a correspondence. 



S 



4.1 Construction of A and bisection algorithm 

A direct consequence of the theorem is the following algorithm to construct A, which mimics the 
one used to construct £: let uj = oj\ . . .uj n -\l be a binary word of length n (ending with 1) and 
d = O.oj the corresponding dyadic rational. Setting uj* = uj\ . . . cj n _iO, we see that d has exactly 
two representations: 

d = 0.ujO and d = 0. uj*1 

where, if n is any finite binary word, rj G {0, 1} N denotes the infinite sequence obtained by repeating 
r\ indefinitely. Now define the dyadic interval associated to d as the open interval Jd = (r~,r + ) 
whose endpoints are the (non-dyadic) rationals^] 

r~ = O.lo* and r + = O.ujuj, 
Example. Take d = 13/16 = ip(2/5). Since the binary expansions of d are 0.11010 and 0.1100T, 



4.2 



we get r~ = 0.1100 = 4/5 and r+ = 0.11010010 = 14/17. (See also the example in Section 
Proposition 4.1. We have 

A=[0,1)\ (J J d (11) 

de(o,i]nQ 2 

with Q 2 := {k/2 s : s G N, k G Z}. 

Proof. From Theorem |1.1| and Equation ^ . □ 

Note that in the previous construction there is substantial overlapping among the Jd- It is possible, 
however, to directly produce all connected components of the complement of A by a bisection 
algorithm: namely, one can produce a new component in the gap between two previously computed 
components. 

Definition 4.2. Let I = [a, b] C [0, 1] be a closed interval, a^b, and let the binary expansion^ of 
a and b be 

a = 0.a\a2 ••• b = O.6162 • • • 
Let n := min{i : a» 7^ b{\. The binary pseudocenter of I is the dyadic rational 

d{I) :=0.oi...a n _il (12) 
By translating one gets the following description of A: 
Proposition 4.3. Let Fq = [0, 1], and for any n define recursively 

Q n := {/ I / is a connected component of F n , L not an isolated point } 

Fn+l ■= F n \ [J J d (l) 

Then 

a= p| f h 

Proof. From Proposition 2.11 of |CT] and Theorem |1.1| □ 

2 Since the correspondence is orientation reversing, the string defining r~ corresponds to the string A + which 
defines the upper endpoint of the interval in [CTj . Section 2.2. 

3 If a or b are dyadic numbers, take the expansion of a which terminates with infinitely many zeros and the 
expansion of b with infinitely many ones 
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4.2 Period doubling 

Period doubling bifurcation is a well known phenomenon which occurs with a universal structure 
in families of unimodal maps. Period doubling cascades are in one-to-one correspondence with 
periodic windows consisting of isolated points in the set A \ {0}, as shown in jAClj . |IP| . We will 
see how similar windows of isolated points occur in £, and in both cases the combinatorial pattern 
can be understood in terms of the Thue- Morse sequence t := (t n )^^ = (0,1,1,0,1,0,0, 1,1,...). 
Recall that t can be defined in terms of the power series expansion 

Y[(l - z 2 ") = l- z -z 2 + z 3 + --- = J2(-l) tn z n (13) 

fc>0 n>0 

Definition 4.4. Define the map A on the set of finite binary words as 

A (7?) := 77 1 7) (14) 
Moreover, for any finite binary word r\, let tq(t]) be the rational number 0.r]0, and let for any j 

T j(v) ■■= t (A j (rj)) 7-00(77) = lim Tjiv) (15) 
The periodic window associated to 77 is 

W v = [t (?/), 7-00(77)) 

Let 1 £ A be a point with periodic binary expansion, and let x = O.77O with 77O minimal period. 
Then the periodic window associated to 7/ contains exactly the sequence of isolated points just 
constructed: 

i>o 

Let us call generating number of the window W v the pseudocenter do(n) of the first interval 
[To(r?),Ti(?7)]. 

Examples. 

If 77 = e, the empty word, we obtain the sequence {Tj(e)}j>o = {0, |, |, j|, . . . } which converges 
to 

t~oo(c) — 0.824908..., and the generating number of W e is do(e) — 1/2. Note that the binary 
expansion of Too(e) = 0.11010011001011... is the Thue-Morse sequence. 

Taking r/ = (11) we get the sequence {Tj(ll)}j>o = {f>|;ff>---} wm ch converges to the irrational 
number Too(H) = 0.892360... and belongs to the next-to-largest window W\\ in A. Its generating 
number is do(ll) = 7/8. 

Windows of isolated points are present also in the set £ and can be generated with a similar 
procedure, introduced in [CMPT] and |CTj without connection to bifurcations of unimodal maps. 
Indeed, let I r be a maximal quadratic interval, and let r = [0; So] = [0; Si] be the continued fraction 
expansions of its pseudocenter, in such a way that So has even length and Si has odd length. Let 
us define So := So, Si := Si and E n +i is generated from S n via the substitutions Si — > Si So, 
So — > Si Si. If we denote 

a n (r) := [0; S n ] a^r) := lim a„(r) 

n— >oo 

it follows immediately from |CT| . Proposition 3.9 that 
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Lemma 4.5. For each I r maximal quadratic interval of pseudocenter r , 

£ n (aoo(r), ao(r)] = (J a n {r) 



n>0 



Limit points aoo( r ) of these cascades have an aperiodic continued fraction expansion with a 
common combinatorial structure which is again related to the Thue-Morse sequence t = (i n )o° : 
indeed, for any r 

«oo(r) = [0; S kl ,Sk 2 ,Sk 3 ,- ..] 

where the sequence K := (ki, fe, ...) = (1,0,1,1...) is such that t(K) = X)j=o For 
instance, if r = \, t herfl a« ,(§) = [0 ; 2, 1, 1, 2, 2, 2, 1, 1, 2, . . . ] ^ 0.38674997... is the accumulation 
point considered in [CMPTj and |KSSj . 

We end this section with a precise characterization of the points {tj(^)}j>o in A and of their 
accumulation point Too (77), which we prove to be a transcendental number. 

Lemma 4.6. Zei n = rji . . . r] p be a finite binary word of length p > suc/i f/mi To(ry) € A and 
define recursively the integers: 

• ^0 := ELi ^2f +1 - 1 

• P i+ i := (P i + 2)(2 2J (f+ 1 )-l) 
. Qj := 2 2J (p +1 ) - 1 

Then Tj(n) £ A for every j, and Tj(j]) = Pj/Qj. Moreover, the binary pseudocenters are given by 



Proof. We need the following elementary identity 



o. ^...^ = ^ (i6) 

x * 2 s — 1 * 2* 

i=l 



for any binary word ui\ ...ui s . Pq is the integer whose binary expansion is 77O, hence from (16) 
To = Po/Qo- Moreover, if P is the integer whose binary expansion is £0, with £ a binary string, 
then the number whose binary expansion is £1£0 is P' = (P + 2)(2^ +1 - 1). By the recursive 
formulas (14) and (15) and identity (16) one has tj = Pj/Qj for all j > 1. Finally, the binary 

pseudocenter of the interval [0.£0, 0.£1£0] is d = 0.£1 = (P + 1)2 Cl^l"^ 1 ) hence the lemma follows 
by taking P = Pj. □ 



Furthermore, let us consider the generating function (cfr. Equation (13)) 

E(z):=H(l-z 2k ) (17) 

This function satisfies the functional equation H(z) = (1 — z)'E(z 2 ), from which we see that all its 
zeroes lie on the unit circle and are actually dense there (see also [Is]). 



4 Notice that the sequence of partial quotients of a ao (^) coincides with the first difference sequence of 



c =13457911 12 13 ■•■ defined as n G c <S=> 2n 4 c. 
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Proposition 4.7. Let 77 = 771 . . . rj p be a finite binary word of length p > such that 70(77) G A. 
Then 



7-00(77) = l-(l-dofo)) 2 
/n particular, T^rj) is transcendental. 



1 



2 p+i 



(18) 



Proof. Setting 7^(77) = Pj/Qj and using Lemma 4.6 we get for j > 1 



Q, - = (2 



(p+l)2J- 1 



i)(Q 



i-i - r i 



Pi-x) 



(Qo-Po)]!^ 2 *- 1 ) 



fc=0 



Therefore 



r iW - 1 " Wo - Po) 2(p+1)2J - _ 1 1 ^7 



2 -2 i (p+l)) 



2 _2J(p+l) 



(l-d0(77)) 



where we used the relation Q - Po = 2 P+1 (1 - do (77)) given by Lemma 4.6 The claimed identity 
follows by taking the limit. Finally, Mahler proved in ([M]) page 363) that 3(a) is transcendental 
for every algebraic number a, with < \a\ < 1. □ 



Examples of Too (77) are 

1/1 

Too(e) = 1 - - i 



Too(H) = 1 



1 _ /l 



and Too (1101) = 1 



5 „ / 1 



32 V 32 



A corresponding transcendence result can be given for £. Namely, accumulation points of period 
doubling cascades a^r) have continued fraction expansion given by the substitution rule explained 
above, hence they are transcendental by [ABJ. 

It has been widely conjectured that all numbers with non-eventually periodic bounded partial 
quotients are transcendental. Let us point out that such a statement is equivalent to the transcen- 
dence of all non-quadratic irrationals in E. Similarly, one might ask whether all irrational points 
in A are transcendental as well. 



4.3 The topology of bifurcation sets 



By using the bisection algorithm and Proposition 4.7 it is immediate to determine the topology of 
bifurcation sets, namely 

Proposition 4.8. Let Qi be the set of rational numbers with odd denominator. Then 

1. The set A n Qi is dense in A. 

2. The derived set A' (i.e. A minus its isolated points) is a Cantor set. As a corollary, A' and 
A have cardinality 2 U . 

3. Let C be the usual Cantor middle-third set, let VLk = (ctk)fik)> k > 1 be the connected compo- 
nents of the complement [0, 00) \ C. For each k, pick a countable set of points Pk Q fife which 
accumulate only on a^. Then A is homeomorphic to C U Ufc>i Pk- 
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4- Any open neighbourhood of an element r G A\Qi contains a subset of A which is homeomor- 
phic to A itself (fractal structure). 

Note that, even though the proposition was stated for A, the same result holds for £ since they are 
homeomorphic. The only change is that you have to replace Qi by the set of quadratic irrationals 
with purely periodic continued fraction expansion. 



Proof. 1. A is closed with no interior by Proposition 4.1 hence the endpoints of the connected 
components of its complement are dense in it. The connected components of the complement 
of A are intervals of type Jd, so their endpoints have purely periodic binary expansion, hence 
they belong to Qi. 

2. By virtue of the homeomorphism cp of Theorem |1.1[ it is equivalent to prove the same result 
for £. Let V denote the points of £ which are periodic under G 

V := {A G £ : G ko (X) = A for some k G N} 
and let us define the set of primitive elements as 

Vo := {A 6 V : A has even minimal period}. 
Lemma 4.9. The set of isolated points of £ is precisely V\Vq. 

Proof. If a is not a limit point of £ then it is the separating element between two adjacent 
maximal quadratic intervals Jo and J\ , so a is the left endpoint of the rightmost interval and 
has odd minimal period (see |CT| ). □ 

Moreover, if A E V, we denote with. cXoc (A) the limit point of the period doubling cascade 
generated^ by A and let W\ := (aoo(A), A) be the corresponding period doubling window. We 



also set Voo '■= {«oo(A) : A G V}, hence by Lemma 4.9 



£' = [0,1)\ (J («oo(A),A). 

Now, no two intervals of the form (aoo(A),A) with A G Vq are adjacent to each other, since 
the right endpoint A is always a quadratic irrational while aoo(A) is always transcendental 



(see Proposition 4.7), therefore £' has no isolated points and it is a Cantor set. 



3. Every closed subset of the interval with no interior and no isolated points is homeomorphic 
to the usual Cantor middle-third set via a homeomorphism of the ambient interval. Such an 
extension can be chosen so that it maps any period doubling cascade to some 

4. Immediate from 3. 

□ 



5 So for instance a 00(g) = 0.38674997... ( |CMPT| . page 24). 
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Let us remark that similar results have been obtained for the set of univoque numbers in |EH Jj . 
[KL2], and by Allouche and Cosnard for a similar number set (PQ, [AClj . [AC2]). For details on 
the relations between these sets see Section [5.21 

Moreover, £ shares some features with the Markov spectrum, since the Gauss map is the first- 
return map on a section of geodesic flow on the modular surface. Even though the two sets are not 
homeomorphic, property 2 also holds for the Markov spectrum [MoJ. 



4.4 Hausdorff dimension 



In |CT] it is proved that dimnS = 1 by providing estimates on the Hausdorff dimension of its 
segments. More precisely, for every K G N+ we consider 

£ K := £ n [1/(K + 1), 1/K], B K := {x = [0; a h a 2 , a 3 , ...] : a* < K Vt G N}, V K (x) := 

K + x 

and it is easy to check that 

*k{B K -i) C£ K C B k , (19) 

therefore]^] diniH£ = limx-^oo dimn^K = ^im-K-Hx dimnBK = 1. 

We can use our dictionary to obtain analogous results in the linear setting, where in fact one can 
explicitely compute the Hausdorff dimension of segments of A. Let Ak '■= An [1 — 2~ , 1 — 2" 
$k(x) '■= x ^ 1 - 2~ K x and 

Ck '■= {x G [1/2, 1] : x does not contain sequences of K + 1 equal digits}; 

using the correspondence of Theorem |1.1[ the inclusions (19) become $k(Ck-i) C Ak C Ck, and 
thus 

dimnCK-i < dimnAK < dimnCK- (20) 

The set is self-similar, therefore its Hausdorff dimension can be computed by standard tech- 
niques (see [F], Theorem 9.3). More precisely, if a^(n) is the number of binary sequences of n 
digits whose first digit is 1 and do not contain K + 1 consecutive equal digits, one has the following 
linear recurrence: Q 

a K (n + K) = a K {n + K - 1) + ... + a K (n + 1) + a K (n) (21) 

which implies 

Proposition 4.10. For any fixed integer K > 2 the Hausdorff dimension ofCx is log 2 (Afc), where 
Xk is the only positive real root of the characteristic polynomial 

P K (t) := t K -{t K ~ l + ... + t+l) 

As a consequence, a simple estimate on the unique positive root of Pk yields 

dimfl- A = lim dimnCx = 1 

K— !-+oo 



6 In fact, the asymptotics of [He] yields the estimate dimu£-K = 1 — 6/iv 2 K + o(l/K). 

7 Sequences satisfying this relation are known as multinacci sequences, being a generalization of the usual Fibonacci 
sequence; the positive roots of their characteristic polynomials are Pisot numbers. 
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5 Kneading sequences, complex dynamics and univoque numbers 



The goal of this section is to explore the relations between the bifurcation sets we described before 
and other well-known sets for which a combinatorial description can be given, namely the real slice 
of the Mandelbrot set and the set of univoque numbers. 

5.1 Relation to the Mandelbrot set 

Let us recall that the Mandelbrot set encodes the dynamical properties of the family of quadratic 
polynomials 

p c (z) := z 2 + c c G C 

namely it can be defined as 

M := {c G C | Pc(°) does not tend to 00} 

A combinatorial model of the Mandelbrot set can be constructed by using the theory of invariant 
quadratic laminations, as developed by Thurston [ThJ. 

Consider the unit disk in the complex plane with the action of the doubling map 

/(*) = * 2 

on the boundary. A leaf £ is a simple curve embedded in the interior of the circle joining two 
points on the boundary. It is usually represented as a geodesic for the hyperbolic metric in the 
Poincare disk model. One can extend the action of / to the whole lamination: the image of a leaf 
C is, by definition, the leaf which connects the images of the endpoints of C We define the length 
£(C) of a leaf to be the (euclidean) length of the shortest arc of circle delimited by its endpoints, 
measured on the boundary circle and normalized in such a way that the whole circumference has 
length 1 (hence, for any leaf C, < £(C) < \). A leaf is real if it is invariant with respect to 
complex conjugation. The doubling map / preserves the set of real leaves. 

A lamination is a closed subset of the disk which is a union of leaves, and such that any two 
leaves in the lamination can intersect only on the boundary of the disk. A quadratic invariant 
lamination is a lamination whose set of leaves is completely invariant with respect to the action of 
the doubling map (we refer to |Thj . page 66 for a complete definition of this invariance). Given 
an invariant quadratic lamination, its longest leaves are called its major leaves, and their image 
the minor leaf. There can be at most 2 major leaves, and they both have the same image. In 
the well-known theory of Douady, Hubbard and Thurston, the leaves of an invariant quadratic 
lamination define an equivalence relation on the boundary circle and the quotient space is a model 
for the Julia set of a certain quadratic polynomial. 

The quadratic minor lamination QML is the set of all minor leaves of any quadratic invariant 
lamination 

QML := {£ I 3 an invariant quadratic lamination Lq s.t. C is its minor leaf} 

By ([ThJ, Theorem II. 6. 8), different minor leaves do not cross, and QML is indeed a lamination. 
Similarly to the Julia set case, the space obtained by identifying points on the disk which belong 
to the same leaf of QML is a model of the classical Mandelbrot set. 
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We will consider the set RQML C QML of real leaves inside the quadratic minor lamination: 
under the previous correspondence, they correspond to real points on the boundary of the Man- 
delbrot set. Moreover, let 1Z = RQML n S* 1 be the set of endpoints of all leaves in RQML. Since 
RMQL is symmetric, it is enough to consider the "upper half" TZn [0, g]. 




(a) The invariant quadratic lamination for the pa- (b) The quadratic minor lamination QML is the col- 

rameter | G A. In red the minor leaf joining ( jjj, j^), lection of all minor leaves. Real leaves are drawn in 

in blue the major leaves. The corresponding param- red, and the set 1Z is the set of all endpoints of red 

eter in £ through the dictionary of SectionWlis 3 ~ 2 x/5 , leaves. A thicker red line represents the minor leaf 

the square of the golden mean. corresponding to the lamination to the left. 



Proposition 5.1. The map x i— > 2x maps TZn [0, \] bijectively onto A. Hence 



K 'K) U M A 



Proposition 5.1, together with the fact that dimnA- = 1, yields a new proof of the main result 
in Za 



Corollary 5.2. The Hausdorff dimension oflZ is 1. 

Because of the fact that the set of rays which possibly do not land has zero capacity and by 
Makarov's dimension theorem [Ma] , we proved 



Corollary 5.3. The intersection of the boundary of the Mandelbrot set with the real line has 
Hausdorff dimension 1 . 



Proof of Proposition 5. 1 . The first ingredient is the fact that the action induced by the doubling 
map on the lengths of leaves is essentially given by the tent map, namely for any leaf C 

2£(f(C)) = T(2£(C)) (22) 
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Let x G TZ fl [0, |] be the endpoint of a minor leaf, let m be the minor leaf and M be one of 
the corresponding major leaves. By maximality of M and invariance of the quadratic lamination, 
£(M) > £{f k {M)) for all k G N, hence by (|2j) 



2£(M) > T k {2£{M)) VA; G N (23) 

From this and the fact that x = 2/3 is a fixed point for T it follows that the major leaf is longer 
than |, hence 

£(m) = -T(2f(M)) = 1 - 2£{M) 
Moreover, by symmetry of the minor leaf with respect to complex conjugation, 

£{m) = l-2x 



hence x = £(M) and by Q 

2x > T k (2x) Vk G N 

hence 2x belongs to A. Conversely, if y = 2x G A, one can construct the real leaf m which joins x 
with 1 — x. The fact that m is a minor leaf follows from the criterion ( |Thj . page 91): a leaf m is 
the minor leaf of an invariant quadratic lamination if and only if 

(i) all forward images have disjoint interior; 

(ii) no forward image is shorter than m; 

(hi) if m is non-degenerate, then m and all leaves on the forward orbit can intersect the (at most 
2) preimage leaves of m of length at least g only on the boundary. 



Indeed, (i) follows since images of real leaves are real; (ii) follows from (23), and (hi) follows 
from the fact that the preimages of a real leaf of length less than | are also real. □ 

5.2 An alternative characterization and univoque numbers 

Definition 5.4. A real number q G (1,2) is called univoque ifl admits a unique expansion in base 
q, i.e. if there exists a unique sequence {ck} G {0, 1} N such that 

Y^c k q- k = l (24) 

k>l 

The set IA of all univoque numbers is called the univoque set. 

In the 90's, }EHJ] discovered that there are infinitely many univoque numbers. In fact [EJK 
showed that, given any binary sequence (cfe)22 =1 , then equation (24) defines a univoque number if 
and only if the following condition is met 

(c fc+ i,c fc+ 2,Cfc+3, •••) < (ci,c 2 ,c 3 ,...) ifcfc = VA; > 1 (25) 
(c k+ i,c k+2 ,c k+3 , ...) > (cx,c 2 ,c 3 , ...) ifc fc = l 

where the minor and major signs refer to the lexicographical order and c = 1 — c as usual. 
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Binary sequences which satisfy equation (25) are called admissible. It is not difficult to establish 



(by means of the so called greedy algorithm) that equation ( 24 ) defines an order preserving bijection 
between IA and the set of all admissible sequences. 

The properties of IA were studied by various authors, leading in particular to the determination 
(in |KLlj ) of the smallest element, q± = 1.78723..., which is not isolated in IA and whose corre- 
sponding admissible sequence is the shifted Thue-Morse sequence {t n ) n >\. The geometric structure 
of IA is specified by the following properties: 

• the set IA has 2 U elements; 

• the set IA has zero Lebesgue measure; 

• the set IA is of first category; 

• the set IA has Hausdorff dimension 1. 

For an account of these and other properties of the set IA see the recent paper [KL2J and references 
therein. 

In a series of papers (see [7Q, |AClj . |AC2] ). J.-P. Allouche and M. Cosnard introduced the 
number set 

r := {x G [0, 1] : l-x< {2 k x} < x, Vk G N} (26) 

and recognized ( |AC2j . Proposition 1) admissible binary sequences to be in a one-to-one correspon- 
dence with the elements of T which have non-periodic expansions. It is not difficult to realize that 
the set r essentially coincides with A, namely 

Lemma 5.5. 

A\{0} = T (27) 



Proof. Let x = O.U1U2 ■ ■ ■ be a binary expansion of x G [0, 1]. From ([2]) it follows that 
T k (x) = 



0.u: k+1 uj k+ 2 ... if Ya=i = (mod 2) 
• • • if EiU UJ i = 1 ( mod 2 ) 



Define 



S k (x) := max{{2 fc x}, 1 - {2 k x}} 



O.ujk+iiOk+2 • • • if ^fc+i = 1 

O.U k+1 LJ k+2 ... if UJ k+ i = 

then r = {i£ [0, 1] : S k (x) < x, Vk G N} and since T k (x) < S k {x) we have T C A. 
Let now x G A \ {0}, and let / = {ifcjfceN = {i G N : T\x) = 0.1 ... } (0 G / since x / 0). Given 
/ G N, either I £ I or I £ I. If I G J, then Si(x) = T l (x) < x. If instead / ^ /, then there exists a 
unique k s.t. i k < I < i k +i. Thus 

T ik (x) = 0. 1^.0... , T l {x) = 0.0^01... 

ik+l—ik ik+\ — l 

and hence 

Si{x) = 0. 1^10. . . < 0. L^IO... = T ik (x) < x 

□ 
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Therefore A can be given the following arithmetic interpretation (see also |AC2j ): 

Proposition 5.6. There is a one-to-one correspondence between U and A\Qx- More specifically, 
a number r = X^>i c k 2 _fc with non-periodic binary expansion belongs to A if and only if 1 
Sfc>i c k a ~ k f or some 1 < q < 2 and the expansion is unique. 



Proof. It is an easy consequence of condition (25) that a sequence (c/A is admissible iff 



(ci,c 2 ,c 3 ,...) < (c k+ l,C k+ 2,Ck+3, ■■■) < (ci,c 2 ,c 3 , ...) \fk>l 

In other words, (c/c) is admissible iff the real value r = X/fc>i c k is a nonperiodic element of T; 
thus our claim follows from Lemma 15.51 □ 



It has been shown in [AC3 ] that q\ £ Li is transcendental. The argument used there can be 
adapted to show that the image in U of all the numbers Too (n) dealt with in Proposition |4.7| are 
transcendental as well. 

Proposition 5.7. Let ^(r/) £ A be an accumulation point of a period doubling cascade, and let 
q G [1,2] be the corresponding univoque number. Then q is transcendental. 

Proof. Let rj be a finite binary word of length p, and write Too(^) = Sfc>i c fc2~ fc and 1 — do(^) = 
Sf=i a «2 - *, with Ck,ai G {0,1}. Define the formal power series A(z) := ^ fc>1 (l — c k )z k and 



:= (aiz + ■■■ + a p z p )E(z p+1 ). By Equation @ = hence A(*) = -B(z) as formal 



power series, since 1—^(77) has a unique binary expansion being transcendental. As a consequence, 
A(-) = and since X)jt>i c fe9 _fc = 1; ^(«) = ~ !• ^ 9 i s algebraic, then A(^) is algebraic, 

hence also and ^(^py) are algebraic, contradicting Mahler's criterion ([Ml) page 363). □ 

From our construction it follows that rational elements of A correspond to algebraic elements in 



U (and are in fact dense in U jdVj ) : rephrasing the question posed at the end of Section 4.2 
might ask whether these are the only non-transcendental elements of U. 



one 



Remark 5.8. Other geometric properties of the set U also follow easily from our correspondence. 
Let u : A \ Qi — > [1, 2] be the map which associates to r = ^^>i 2~ fc the unique value q which 



satisfies the relation (24). Lt is not difficult to check that u is locally Holder- continuous and extends 
to a homeomorphism between A' and the closure U of U in [1,2], hence U is also a Cantor set. 
Moreover, one can prove that the Holder exponent of u restricted to the segments Ak ( considered 



in Section 4-4) greater than the Hausdorff dimension of Ax, thus dimn u{Ak) < 1 for all K > 1 



and IA has zero Lebesgue measure. 
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